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Abstract 

As the third stage of the project multi-indexed orthogonal polynomials, we present, 
in the framework of 'discrete quantum mechanics' with pure imaginary shifts in one di- 
mension, the multi-indexed Wilson and Askey- Wilson polynomials. They are obtained 
from the original Wilson and Askey- Wilson polynomials by multiple application of the 
discrete analogue of the Darboux transformations or the Crum-Krein-Adler deletion of 
'virtual state solutions' of type I and II, in a similar way to the multi-indexed Laguerre, 
Jacobi and (g-)Racah polynomials reported earlier. 

1 Introduction 

This is a third report of the project multi-indexed orthogonal polynomials. Following the 
examples of multi-indexed Laguerre and Jacobi polynomials [TJ, multi-indexed (g-)Racah 
polynomials [2], we present multi-indexed Wilson and Askey- Wilson polynomials constructed 
in the framework of discrete quantum mechanics with pure imaginary shifts [3j. It is well- 
known that the original Wilson and Askey- Wilson polynomials are the most generic members 
of the Askey scheme of hypergeometric orthogonal polynomials [HEJEIE]- These new multi- 
indexed orthogonal polynomials are specified by a set of indices T> = {di, . . . , 6m} consisting 
of distinct natural numbers dj G N, on top of n, which counts the nodes as in the ordi- 
nary orthogonal polynomials. The simplest examples, T> = {£}, £ > 1, {P^ n (x)} are also 
called exceptional orthogonal polynomials [8]-[29|. They are obtained as the main part of 



the eigenfunctions (vectors) of various exactly solvable Schrodinger equations in one dimen- 
sional quantum mechanics and their 'discrete' generalisations, in which the corresponding 
Schrodinger equations are second order difference equations [3j [30j [31] . They form a com- 
plete set of orthogonal polynomials, although they start at a certain positive degree (£ > 1) 
rather than a degree zero constant term. The latter situation is essential for avoiding the 
constraints of Bochner's theorem [32] . We strongly believe that these new orthogonal poly- 
nomials will find plenty of novel applications in various branches of science and technology 
in the good old tradition of orthogonal polynomials. 

The basic logic for constructing multi-indexed orthogonal polynomials is essentially the 
same for the ordinary Schrodinger equations, i.e. those for the Laguerre and Jacobi polyno- 
mials and for the difference Schrodinger equations with real as well as pure imaginary shifts, 
i. e. the (g-)Racah polynomials and the Wilson and Askey- Wilson polynomials, etc. The main 
ingredients are the factorised Hamiltonians, the Crum-Krein-Adler formulas [331 EH [35] for 
deletion of eigenstates, that is the multiple Darboux transformations [36] and the virtual 
state solutions pQ which are generated by twisting the discrete symmetries of the original 
Hamiltonians. Most of these methods for discrete Schrodinger equations had been developed 
[30| [26l [3l [311 E3 EHl [39] and they were used for the exceptional Wilson and Askey- Wilson 
polynomials [121 EH] ■ 

It is important to stress that the factorised Hamiltonians in the discrete quantum mechan- 
ics, that is, those governing the (g-)Racah, Wilson, Askey-Wilson polynomials etc possess 
certain discrete symmetry. They lead to virtual Hamiltonians which are linearly connected 
with the original Hamiltonian. (See fl2.14|) - fl2.16l) of the present paper and (2.18)-(2.22), 
(2.59)-(2.63) of [2].) In the ordinary quantum mechanics of the radial oscillator potential 
(x 2 + ^11, for the La guerre polynomials) and the Poschl- Teller potential ( . $ + h( * h 2 - , 
for the Jacobi polynomials), the discrete symmetry is well known. For the former, g — > 1 — g 
and/or x — > ix and for the latter g — > 1 — g and/or h — y 1 — h. To the best of our knowl- 
edge, those discrete symmetries for the (g-)Racah, Wilson and Askey-Wilson systems do 
not seem to be widely recognised, since they are not easily identifiable in the polynomial 
equations. The virtual state solutions belong to the virtual Hamiltonians. The actual con- 
tents of virtual state solutions depend on the types of the Schrodinger equations. For the 
ordinary Schrodinger equations with a second order differential operator, the virtual state 
solutions satisfy the Schrodinger equation. But they do not belong to the Hilbert space of 



2 



square integrable solutions, due to the twisted boundary condition on either one of the two 
boundaries, to be called the type I or II. For the multi-indexed (g-)Racah polynomials in 
the discrete quantum mechanics with real shifts, the virtual state 'solutions' fail to satisfy 
the Schrodinger equation at either one of the two boundary points [2]. They are called vir- 
tual state vectors of type I or II. In the present case of discrete quantum mechanics with 
pure imaginary shifts, the virtual state solutions satisfy the difference Schrodinger equation. 
But they do not belong to the Hilbert space of eigenfunctions either by the lack of square 
integrability or by the presence of singularities in certain rectangular domain. (See more 
detailed discussion in section 12.21 ) In other words, the analyticity requirements supersede 
the boundary conditions which used to classify the virtual state solutions for the Laguerre, 
Jacobi and (g-)Racah cases. In section three, we will introduce two types of twistings or the 
discrete symmetry transformations and the corresponding virtual Hamiltonians and virtual 
state solutions. They are of the same structure but adopting different sets of parameters. 
We will call them of type I and II as in the other cases but they are not related to boundary 
conditions. In all these cases, the features disqualifying them to become the eigenfunctions 
are carried by the so called "virtual groundstate" functions 4>o(x), (I3.28p . The polynomial 
part of the virtual state solutions, to be denoted by £ v (??) (13.281) . are the genuine solutions 
of equations determining the eigenpolynomials, but with twisted parameters. It is the vir- 
tual state polynomials {£ v (77)}, not the virtual groundstate (fio(x), that play the main role 
in the construction of multi-indexed and exceptional [T2j [19] polynomials and the set of 
their degrees {di, . . . , dm} constitutes the multi-index . We focus on the algebraic structure 
of the multi-indexed orthogonal polynomials and their difference equations, which hold for 
any parameter range. We do not pursue the other important aspect of the problem, that 
is the determination of the parameter ranges in which the hermiticity of the multi-indexed 
Hamiltonians and the positivity of the orthogonality weight functions for the multi-indexed 
polynomials are ensured. 

This paper is organised as follows. In section two, the basic logic of virtual states dele- 
tion in discrete quantum mechanics with pure imaginary shifts in general is outlined. Start- 
ing from the general setting of discrete quantum mechanics with pure imaginary shifts in 
§ 12. 14 the analyticity requirements in connection with the hermiticity (self-adjointness) of the 
Hamiltonians are briefly recapitulated in § 12.21 General procedures and formulas of multiple 
virtual states deletion are reviewed in § 12.31 The main logics are essentially the same as 
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those for the multi-indexed Laguerre, Jacobi and (g-)Racah polynomials but the explicit 
formulas look rather different reflecting the specific properties of discrete quantum mechan- 
ics with pure imaginary shifts. After recapitulating the basic properties of the Wilson and 
Askey- Wilson systems in § 13. 1[ the discrete symmetries of the the Wilson and Askey- Wilson 
systems are introduced in § 13.21 The multi-indexed Wilson and Askey- Wilson polynomials 
are constructed explicitly in § 13.31 for the type I and II virtual states deletions. The ana- 
lyticity and hermiticity of the multi-indexed hamiltonians are discussed in some detail in 
section §3.41 The final section is for a summary and comments including the limits to the 
multi-indexed Jacobi and Laguerre polynomials. For simplicity of presentation we relegate 
several technical results to Appendix. 

Throughout this paper we will focus on the algebraic aspects of the theory. To determine 
the exact ranges of validity of various formulas is another problem. 

2 Formulation 

2.1 Original system 

Let us recapitulate the discrete quantum mechanics with pure imaginary shifts developed 
in [3]. The dynamical variables are the real coordinate x (x\ < x < X2) and the conjugate 
momentum p = —id x , which are governed by the following factorised positive semi-definite 
Hamiltonian: 

U = ^/V{^e <p ^V^ + ^V^e-^^V{^ -V{x) -V*{x) = A j A, (2.1) 
A = i(e&y/V^-e-&y/V@j), ^ = ~i(Vv(xj e * P ~ ^/V*{x) e~i p ). (2.2) 

Here the potential function V(x) is an analytic function of x and 7 is a real constant. The *- 
operation on an analytic function f(x) = G C) is defined by f*(x) = 'Yl, n Oj* n x n 1 

in which a* is the complex conjugation of a n . Obviously f**(x) = f(x) and f(x)* = f*(x*). 
If a function satisfies /* = /, then it takes real values on the real line. Since the momentum 
operator appears in exponentiated forms, the Schrodinger equation 

U<P n (x) = £ n <P n (x) {u = 0, 1 , 2, . . .) , (2.3) 

is an analytic difference equation with pure imaginary shifts instead of a differential equation. 
Throughout this paper we consider those systems which have a square-integrable groundstate 
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together with an infinite number of discrete energy levels: = £q < £\ < S 2 < ■ ■ ■ ■ The 
orthogonality relation reads 



X2 

dx <p* n (x)<f) m (x) = h n 5 nm (n,m = 0, 1,2, .. .), < h n < 00. (2.4) 

XI 



The eigenfunctions (f) n (x) can be chosen 'real', = 0„(x), and the groundstate wave- 

function 4>o(x) is determined as the zero mode of the operator A, A4>o(x) = 0, namely, 



yJv*{x-q)Mx - ii) = y/v{x + q)Mx + ^ 2 L )■ (2-5) 
2.2 Analyticity requirements 

The hermiticity of the Hamiltonians of discrete quantum mechanics with pure imaginary 
shifts is more involved than that of the ordinary quantum mechanics [31 [3TJ . Here we review 
the hermiticity of the Hamiltonian (12. ip in general, in a way applicable to those appearing in 
the multi-indexed Wilson and Askey- Wilson systems, e.g. fl2T8|) - fl2T20|) . fl2T26|) -f l2T28|) . Of 
course, the hermiticity of the original Hamiltonians of the Wilson and Askey- Wilson systems 
( jSaD - fllDID is well established [31 ETJ. 

Let us consider the functions of the form f(x) = 4>q(x)1Z(x) , where (po(x) 2 and TZ(x) 
are meromorphic functions and 7Z*(x) = 7Z(x). For two such functions j\ = <PqR-i and 
f 2 = 0o^-2, the condition of the hermiticity {fx^Hf^) = (7~t-fi, $2) becomes 

[ 2 dx(G(x-i%) + G*(x + q)) = I 2 dx(G(x + q) + G*(x-q)), (2.6) 

J X\ J X\ 

where G[x) is defined by 

G{x) = v{x + %i)(po{x + q ) 2 7^i(s + q )n 2 {x - q), 

(=> g*(x) = v(x + q)4> ( x + *i) 2 ^i(x - q)n 2 (x + qj) . (2.7) 

Although the term V(x) + V*(x) in % are canceled out in this calculation, this term V(x) + 
V*(x) should be non-singular for x\ < x < x 2 , and the integral f* 2 dx V \x)(j>o(x) 2 'Ri(x)'R, 2 (x) 
should be finite. By using the residue theorem, the condition (12.61) is rewritten as 
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dx(G{x 2 + ix) — G*(x 2 — ix)) — / dx(G{xi + ix) — G*(x\ — ix)) 



2n f] E ^x a {G(x)-G*(x)), (2. 

xo:pole in D 1 



5 



where the residue of the function G(x) — G*(x) is taken at the poles in the rectangular 
domain D 1 : 

D 1 = {xeC\x 1 < Rex < x 2 , |Imx| < ||7|}. (2.9) 

In our previous work on the exceptional Wilson and Askey- Wilson [T2l [T9] , we required 
that G and G* have no poles in the rectangular domain D 7 , as a sufficient condition for 
the hermiticity of the Hamiltonian. This was too strong a requirement. In later examples, 
TZ(x) = TZ[T](x)) is a rational function of i](x) and V(x)(po(x) 2 has the form ~ (^0o) original x 
(rational function of r](x)). In the original Wilson and Askey Wilson theory, (V0q) original 
part has no poles in D 7 , (13.141) . In the deformed theory in general, however, (^0o) O rigmai has 
shifted parameters and we continue to require that this part has no poles in the rectangular 
domain D 1 . We also remark that even if the (rational function of 77(x))-part has poles in the 
rectangular domain D 1 , there is a possibility that the sum of the residues vanish. 

2.3 Deletion of virtual states 

In [3H] we have presented the Crum-Adler scheme, i.e. the deletion of M eigenstates. In 
that case the index set of the deleted eigenstates V = f {di,<i 2 , . . . ,o?m} (dj G Z>o) should 
satisfy the condition U-=i( m ~ dj) > (Vm e Z> ), eq. (2.8) in [38]. We now apply the 
Crum-Adler scheme to virtual states instead of eigenstates. The above condition eq. (2.8) in 
[38] is no longer necessary. 

The Casorati determinant of a set of n functions {fj{x)} is defined by 

W 7 [/ l5 . . . , f n }(x) ^ det (Mxf)) , xf ] d = f x + i(*±l - j)7, (2.10) 

\ / l<j,k<n 

(for n = 0, we set W 7 [-](x) = 1), which satisfies identities 

w 7 [a, . . . , f n \*(x) = w 7 [/r, . . . , r n }(x), (2.ii) 

n 

W^g^gh, ...,gf n } = U^f) ' W t[/i, h, ■ ■ ■ > /«](^). ( 2 - 12 ) 

i=i 

W 7 [W 7 [/ a , / 2 , . . . , f n , g], W 7 [/ ls / 2 , . . . , /„, ft] ] (ar) 

= W 7 [/i, / 2 , • • • , fn)(x) W 7 [/i, / a , . . . , /„, ^, ft](z) (n > 0). (2.13) 
Let us assume the existence of an analytic function V'(x) of x satisfying 
V(x)V*(x - i 7 ) = aY'(i)f(i - i 7 ), a > 0, 
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V(x) +V*(x) = a(V(x) + V'*(x)) -a', a'<0, (2.14) 
where a and a' are constants. Then we obtain a linear relation between two Hamiltonians: 



H = olH! + a', 



(2.15) 
(2.16) 



Since % is positive semi-definite, %' is obviously positive definite and it has no zero-mode. 
Let us also assume the existence of virtual state wavef unctions <f) v (x) (v G V), which are 
'polynomial solutions' of degree v of the Schrodinger equation 



H(f> v (x) = £ v <p v (x) or H'(f) v (x) = £' v <f) v (x), £ v ^a£' v + a', = 4> v (x), (2.17) 



but they, including the zeromode 0o> do not belong to the Hilbert space of H. Here V is 
the index set of the virtual state wavef unctions. We require that £ v < and some analytic 
properties of (f) v (x), which are explicitly presented in §[3J 

We have developed the method of virtual states deletion for ordinary quantum mechanics 
in pQ and for discrete quantum mechanics with real shifts in [2J. Algebraic aspects of this 
method are the same and can be applied to discrete quantum mechanics with pure imaginary 
shifts. The procedure is as follows; (i) rewrite the original Hamiltonian as H = A dl A dl 
(di £ V), (ii) define a new isospectral Hamiltonian % dl = f A dl A d + £ dl , whose eigenfunc- 
tions are given by <f) din (x) = A dl 4> n (x) together with virtual state wavefunctions <fi diw (x) == 
A dl 4> v {x) (v G V\{di}), 'Hd 1 (f>d 1 n{x) = ^(pdmix), 'Hd 1 ^d 1 w{x) = £ v 4> dl v{x) (iii) rewrite this as 
Hdi — A^ dld2 A dld2 + £~d 2 [d 2 G V\{di}), (iv) define the next isospectral Hamiltonian "H<fid 2 
^-d!d 2 A did2 + £ d , 2 , whose eigenfunctions are given by 4> dld2n (x) = A dld2 (j) din (x) together with 
virtual state wavefunctions 4> dld2V (x) = f A dld2 (j) div (x) (v G V\{di,d 2 }), T-L dld2 (f) dld2n (x) = 
Zn&didinix), Wdidi&didivix) = £v<Pd 1 d 2 v(x), (v) by repeating this process, we obtain U dl ...d s 
and its eigenfunctions 4>d 1 ...d s n(x) together with virtual state wavefunctions <j>d 1 ..A a \{%)i ( y i) 
r Hd l ...d a can be written in the standard form, T-i dl ... da = A di dg A dl ... da . If the resulting system 
is well-defined, we obtain the isospectrally deformed systems just as those in refs.pQ and [2j. 

Here we present 'formal' expressions of the deformed systems, which are proved induc- 
tively. The system obtained after s virtual state deletions (s > 1), which are labeled by 
{di, . . . , d s } (dj G V : mutually distinct), is 



7~Ld\...d 3 Ad\...d 3 A ( 



+ £ d . 



(2.18) 
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A dl ... ds =i{e*\IW^x)-e-&,Jv^(x)), 

A... ds = -*(\/v dl ... d3 (x)el p - ^U(*)e-**), (2.19) 
V dl ...d 3 (x) d = f yJv(x-i^)V*(x-i^) 

x W 7 fc 1 ,...,0 da _ 1 ](x + ^) W 7 fe 1 ,...,0 da ](g-H) ^ 2Q ^ 

W 7 [0 dl ,...,0 da _J(x-^) W 7 [0 dl ,...,0 d J(x) 

<?W..d s n(z) = f A dl ...d B <t>d 1 ...d a -in(x) (n = 0, 1,2, . . .), 

4> dl ... ds v(x) = A dl ... da (f) dl ... d3 _ 1 v(x) (v e V\{di, (2.21) 

y. dl ... ds (j> dl ... da n{x) = £n4> dl ... ds n(x) (n = 0, 1,2, . . .), 

'H dl ... ds 4>di...d s v{x) = £v<Pdi...d.v(x) (vG V\{di,...,d s }), (2.22) 

s 

{<f>d 1 ...d s n,<f>d 1 ...d s m) = \\{£n ~ £ dj ) • K5 nm (n, 771 = 0, 1, 2, . . .). (2.23) 

3=1 

Let us remark that the eigenfunctions and the virtual state solutions in all steps are 'real' by 
construction, 0^...^ „(x) = 0^...^ n(^), 0^...^ v( x ) = ( Pd 1 ... ds v(x) and they have Casoratian 
expressions: 

(j) dl ... ds n(x) = A(x)Wj[4) dl , . . . , 4> ds , (j) n ](x), 

) = A(x)W 7 [0 dl ,...,0 ds ,0 v ](:r), (2.24) 

i 



W 7 [0di, • • • , 0dj(^ - »i)W 7 [0di, • • • , 0d(z + i£) 



which are shown by using (I2.13p . 

Writing down (I2.22p and dividing it by 4> dl ... ds n{.x) or 4> dl ... ds and using (12.241) . we 
obtain 

v dl ... ds (x + q) + vi.. da (x - q) - 4 + (2.25) 

= - »§7) ~ ^7) WSd1 ' • • • ' + — ^ An]iX - %1) 

V " " W 7 [0 dl , . . . ,4> ds ]{x - ii) W 7 [0 dl ,... ) ds) n ](a;) 

V W 7 [0 dl ,...,0 d J(x + ^) W 7 [0 dl ,...,0 ds ,0 n ](x) 

and a similar equation for the virtual state solution. 

The deformed Hamiltonian % dl ,., ds can be rewritten in the standard form: 

T-i dl ... da = A dlds A du .. ds , (2.26) 
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A dl ... ds = i {efryjv^ (x) - e-fry/V*...*. (x) ) , 



(2.27) 



V dl ... ds (x) = yJv{x-i^)V*{x-i^) 



W 7 fe, . . . , da ](x + ij) W 7 [0 dl , . . . , da , </> ](g - ij) 
W 7 [0 dl , . . . , 4> da ] (x - %D W 7 [4> dl ,...,<f> ds ,<t> ] (x) 



(2.28) 



in which the A operator annihilates the groundstate, 



A dl ... ds 4> dl ... ds o(x) = 0. 



(2.29) 



The conditions for the equality of (I2.26|) and (12. 18)) are 



V dl ...d s { x ) v L..d s { x - ij) = V dl ... ds (x - i^)Vl___ da (x - q), 

v dl ... da ( x ) + v l...d s ( x ) = v dl ... da ( x + »i) + v£...d s ( x - «i) - £d s 



(2.30) 



The first equation is trivially satisfied and the second equation is a consequence of (12 .25)) 
with n = 0. 

It should be stressed that the above results after s-deletions are independent of the orders 
of deletions {<f) dl ... ds n{x) and 4> dl ... dsV (x) may change sign). 

In order that this deformed system is well-defined i.e. Hamiltonian T-L dl ,„ da is hermitian, 
we have to study the singularities of V dl ,„ ds (x) and 4> dl ... dsn ( x )- We will do this for explicit 
examples in section §3.41 



3 Multi-indexed Wilson and Askey- Wilson Polynomi- 



al this section we apply the method of virtual states deletion to the exactly solvable systems 
whose eigenstates are described by the Wilson (W) and Askey- Wilson (AW) polynomials. 
We delete M virtual states labeled by 



and denote H dl ... dM , <j) dl ... dM n, A dl ... dM , etc. simply by Ut>, <hn, Ax>, etc. 

We follow the notation of [3J. Various quantities depend on a set of parameters A = 



als 



V = {di, d,2, ■ ■ ■ , (1m} (dj G V : mutually distinct), 



(3.1) 



(Al.A: 
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3.1 Original Wilson and Askey- Wilson systems 

Let us consider the Wilson and Askey- Wilson cases. Various parameters are 



•l 1 i i\ 

,2' 2' 2' 1'i 



K 



W : X\ — 0, X2 — oo, 7 = 1, A = (ai, a^, a 3 , a^), <5 
AW: x x = 0, x 2 = 7T, 7 = log 9, g A = (ai,a 2 ,a 3 ,a 4 ), £ = (5,5,5,5), k = g~\ (3.2) 

where g A stands for g( Al > A2 '---) — (g Al ; g A2 ; . . .) and < g < 1. The parameters are restricted 
by 

{a*, a^, 03, a|} = {ai, a 2 , 03, 04} (as a set); W: Reaj > 0, AW: |aj| < 1. (3.3) 
Here are the fundamental data [3J: 



(2ix(2za; + 1)) 1 nj=i(% + 



ar 



: W 



V{x; A) 
77(3;) = 

£»(A) = 

<f) n (x; A) 

P n (x;A) = P„(r ? (x);A) 



cosx : AW ' 
n{n + bi — 1) 



•_]_(%• + ix) 


: W 


)) _1 nj=i(i 


- aj e ix ) : AW ' 


f 2x 


: W 


|^ 2sinx 


:AW ' 


: W 


u def 

Ox = ai + a 2 + a 3 + a 4 


:AW ' 


L def 

O4 — 0x02*33 0-4 , 



X 4 </> £ 



Cn(A) = 

0o (a; A) 



(g- n -l)(l-6 4 g n_1 
^o(x; A)P n (:r; A), 

W n (77(x);ai,a 2 ,a3,a 4 ) :W 
p n {r](x);a 1 ,a 2 ,a 3 ,a i \q) : AW 

(ai + a 2 ) n (oi + a 3 ) n (a 1 + a 4 )„ 

/— n, n + bi — 1, ai + ix, ai — zrr 

X 4.F3 ( 

V ai + a2, ai + a 3 , ai + 04 
ax n (aia2, CI1CI3, aia4 ; q) n 

q~ n , b±q n ~ x ', axe", die""' 

31 

V CL\Qj2i CL\CL 3l O1O4 

c n (A)r](x) n + (lower order terms), 
(_!)»(„ + 6j _ l) n ; W 
2»(6 4 g»-i ; q ) n : AW ' 

(r(2ix)r(-2zx))- 1 nL r(a j + zx)r( aj - %x) 



w 



AW 



: W 



K{\) 



(3.4) 

(3.5) 

(3.6) 
(3.7) 



(3.f 



(3.9) 
(3.10) 

(3.11) 



(e 2i * ; g)oo(e" 2 " ; g)oo E^iK^ ; g)- 1 ^-* ; g)" 1 : AW 
27m! (n + 61 - l) n rii<i<i<4 r (™ + «i + a i) • r ( 2 ™ + biY 1 ■ w 
27r(M n -^g)n(M 2 ";g)oo(g n+ ^g)^ 1 ni< J < J <4(«^^^g)^ 1 : AW ' 

Here W n and p n in (13. 8 j) are the Wilson and the Askey- Wilson polynomials [6] and the sym- 
bols (a) n and (a; g) n are (g-)shifted factorials. We have 0q(x; A) = (po(x; A) and P^(x; A) = 
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P n (x; A). Note that 



{x; X + 6) = <p(x)y/V(x + q; A) o (x + i\\ A), (3.12) 

ip(x) 



V(x;X + d) = k- 1 ^ _ f J r, V(x-q-,X)- (3.13) 



The sinusoidal coordinate r](x) has a special dynamical meaning [30] El ED]- The Heisenberg 
operator solution for r](x) can be expressed in a closed form and its time evolution is a 
sinusoidal motion. The hermiticity of the Hamiltonian is satisfied because the function V(pQ 
has the property 

03. 3 p V(x + q; X)(po(x + q; X) 2 has no poles in the rectangular domain D 7 , (3.14) 

and the function G (\2.7\i (TZ(x) = P n (x] A)) satisfies G{x\ + ix) = G*{x\ — ix), G(x 2 + ix) = 
= G*(x 2 - ix) for W and G(xi + ix) = G*(x\ - ix), G(x 2 + ix) = G*(x 2 - ix) for AW. 
Note that eq. fl3.12l) implies 



V(x + q- A)0 o (x + q; A) 2 = yuv ' ; I ' . (3.15) 

The system is shape invariant [HJ [3], 

^(A)^(A) f = kA{X + 5)U(A + 8) + £i(A), (3.16) 

which is a sufficient condition for exact solvability and it provides the explicit formulas for 
the energy eigenvalues and the eigenf unctions, i.e. the generalised Rodrigues formulas [3]. 
The action of the operators >4(A) and v4.(A)t on the eigenfunctions is 

A(X)<p n (x- X) = / n (A)0„_x (x- X + 6) , A(X)^(p n -i (x; X + S)= 6 n _ 1 (A)0 n (x; A). (3.17) 

The factors of the energy eigenvalue, f n (X) and 5„_i(A), S n (X) = f n (X)b n _i(X) , are given 
by 

/n(A) = | ^_+_ &1 1)(1 1) _ Mn _ 1) [ J w , 6 n _i(A) = | q \ [ J w . (3.18) 

The forward and backward shift operators ^(A) and B(X) are defined by 

J-(A) = o (x; A + S)- 1 o A(X) o o (s; A) = i(p(x)-\e& - e"^), (3.19) 
B(X) = <p a (x;X)- 1 oA(X) ] o0 o (x;A + <5) = -z(F(x; A)e^ p - V*{x; A)e^ p )<^(x), (3.20) 
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and their action on the polynomials is 

F(X)P n (x; A) = / n (A)iVi(z; A + 8), B(X)P n ^(x; X + 8) = b n _ 1 (X)P n (x; A). (3.21) 

The second order difference operator "H(A) acting on the polynomial eigenfunctions is square 
root free. It is defined by 

H{X) = o (x; A)" 1 o H{\) o o (x; A) = B{X)JF{X) 

= V{x; A)(e 7P - 1) + V*{x; A)(e" 7P - 1), (3.22) 
U{X)P n {x; A) = S n (X)P n {x; A). (3.23) 

3.2 Discrete symmetries of the Wilson and Askey- Wilson systems 

Since the potential function V(x; A) (13.41) is invariant under the permutation of (ai, a 2 , a^, a.4), 
the system is symmetric in (01,02,03,04). So are the groundstate wavefunction (f) Q (x; X) 
(I3.10p and the eigenpolynomial P n (x; A) (13. 8ft . 

In the following we restrict the parameters as follows: 

ai,a 2 6l or a 2 = a± ; 03,04 el or a* A = 03. (3.24) 

Let us introduce twist operations t 1 , t n and constants 6,6, 

t : (A) = (1 — Ai, 1 — A 2 , A 3 , A 4 ), 8 = (— |, — 

t n (A) d = f (A 1; A 2 , 1 - A 3 , 1 - A 4 ), 5 H d = f (|, |, -§, -§). (3.25) 
By using these twist operations, we define deformed potential functions V by 

V'\x; A) = V(x; i\X)) , V' l \x; A) = V(x; t n (A)) . (3.26) 

These V's satisfy (f2~T4D with 

i m / 1 :W /i/.x / -(ai + a 2 - l)(a 3 + a 4 ) :W 

" l J \ axasg- 1 : AW ' a [ ' \ -(1 - ai^g- 1 )^ - a 3 a 4 ) : AW ' 

= I 1 1 : IL > «' n ( A ) = I "I? + " 4 " Vi? + : 1, • (3-27) 

\ a 3 a 4 g _1 : AW ' I -(1 - a 3 a 4 g _1 )(l - a a a 2 ) : AW v ; 

(For a > 0, we have aia 2 > and 0304 > for AW case. For a' < 0, we will restrict the 
parameters further as in ( I3.36[) .) We obtain a linear relation between the two Hamiltonians 
(I2.15p . The virtual state wavefunctions satisfying (I2.17P are given by 

ftfoAJ^toMW), A) = f v (x; t : (A)) = 0q(x; A)£*(x; A) (v G V 1 ), 
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= 0v(^A)= f v (x;t n (A)) =ffi(a ? ;A)^(ar;A) (v G V 11 ), 

A) = ^(^(x); A) d = f P v (a;; t n (A)) = P v (^); t n (A)). (3.28) 

The virtual state polynomials £ v (??; A) are polynomials of degree v in rj. They are chosen 
'real,' <p* {x; A) = (j)o(x; A), £,*(x; A) = £ v (x; A) and the virtual energies are £y(\) = £ v (t(A)): 

i/ x x _ / — (a x + a 2 - v - l)(a 3 + a 4 + v) : W 



(3.29) 



1 -(1 - a 1 a 2 q- v - 1 )(l - a 3 a 4 q v ) : AW 

oil /\\ _ f ~( a 3 + a 4 - v - l)(oi + a 2 + v) : W 
v 1 j \ -(1 - a 3 a 4 g^ 1 )(l - aia 2 g v ) : AW ' 

Note that a'(X) = S (X) < and 

W: £j(A) < ai + 02>v + l, £"(A) < <^> a 3 + a 4 >v + l, 
AW : £j(A) < ^ < a x a 2 < g v+1 , £* ! (A) < < a 3 a 4 < g v+1 , (3.30) 

for v > 0. We choose V 1 and V 11 as 

V 1 = {1, 2, . . . , [Ai + A 2 - 1]'}, V 11 = {1, 2, ... , [A 3 + A 4 - 1]'}, (3.31) 

where [x]' denotes the greatest integer not equal or exceeding x. We will not use the label 
states for deletion, see (I3.58p - fl3.59p . 

For later use, we define the following functions (recall Xj in fl2.10p ): 

I, . n def 4>o(x] A) j (M) , def 1/1 ( x j ) ! A ) M , . . 

z/(z;A) = rj(x} ; ;A,M) = — ^ (1 < j < M , 

<feA) v l (x;X+{M -1)6 ) 

ll t , del <Po(x] X) u (M) , b^n def V U (Xj ^ A) . qq\ 

0o(^;A) z/ n (x;A + (M-l)(5 ) 

Their explicit forms are 



A, M) =a I (X+(M- i)tf)-* M -v K HW-V'-V-VW-J) (3.33) 
Yl(a k - M f 1 + ix)j-i(a k - M = 1 - ^) A/ _ J : W 

fc=l,2 

^(M+l^i) J| (^--l^.^^^-^^i,.^^. ;AW > 
fc=l,2 

A, M) = a II (A + (M-l)5 II )- |(M " 1) « |(M - 1)2 ~ a " 1)(M " J ' ) (3.34) 



x 
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FI ( a * ~ ^Y 1 + ix )j~i( a k ~ ^ ~ ix )M-j ■ W 



X 



fc=3,4 



e WM+i-id) ] J (a.g-^e"; ?) i _i(a Jk g— ^e"*; g) M -, : AW 



fc=3,4 



The auxiliary function <fu( x ) [SB] is denned by: 

M-2 

(p M (x) = p{x) [ ^ ] Y\ ((p(x - i\l)y{x + if 7)) 



fc=i 

and (fo(x) — (pi(x) — 1. Here [sc] denotes the greatest integer not exceeding x. 

3.3 Explicit Forms of Multi-indexed Wilson and Askey- Wilson 
polynomials 

We delete M = Mj + M n virtual states V = {d\, . . . , d\ Iv df, . . . , d^}, where d) and d* 1 
are labels of the type I and II virtual states, respectively. We restrict the parameter range 
further as follows: 

W: Recti > |(max{dj} + 1) (i = 1, 2), Re a, > i(max{dj 1 } + 1) (i = 3,4), (3.36) 

AW : |0i| < q h(^Ad}}+l) y = 1)2 ^ | a .| < g I(ma Xj {^}+l) ^ = 3)4 ^ aia2;fl3a4 > 0. 

Then the condition £ V (A) < (v G X>) is satisfied, see f)3.30p . We assume that the parameters 
are so chosen that £ v (x; A) 7^ for x\ < x < x<i- For example, if we take the parameters as 

o m £ R, -KAj--A fc <l (j = 1,2; fc = 3,4), 
o d, a 2 G R, 04 = o 3 , Ai + A 2 + 2v < A 3 + A 4 , 
o a 3 , a 4 el, a 2 = ai, A 3 + A 4 + 2v < Ai + A 2 , 

then £ v (x; A) and £ v (x; A + S) have a definite sign for real x. 

Let us write down <pj) n (12.241) concretely. The Casoratians in f)2.24p are reduced to the 
following determinants, by which we define two polynomials 3d (a:; A) and Pv,n{x', A), to be 
called the denominator polynomial and the multi-indexed orthogonal polynomial, respec- 
tively: 



y(M) f(M) p(M) V?(M) 

°-l Mt "l A'/tt 



(p M {x)E v (x; A) x A, 
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ix, a k - 



M-l 



x Yi k =i, 2 uf=\ 1 (a k 



M-l 



A= ( 



%0C • 



2 lX )3 
M-l 



X 



2 ""J 3 

e ,a k q 2 e ,q) 

M-l ■ M-l 

a k q 2 



: W 



e- ix ;q)j : AW 



(3.37) 



x 2 



\M(M+l) 



^(M+l) 
X 4 



(M+l) y?(A/+l) 



4 1 



y 



(M+l) ^(M+l) 



Mjj 



<^M+l(^)A?,n(^; A) X B, 



B 



M 



IX 



X n fe =i,2 llS K -f + ix,a k -f- ix) 



: W 



Il^ 3A Uf2 1 a- k j q^ +1 \a k q 
X IL=i 2 nS a^qbV+Vfaq-Ze**, a k q-^e^ q)j : AW 



-ix \ 

2 e , a^g 2 e ; q)j 



(3.38) 



where 



=rf(x ( f I ) ;X,M)e v (x^>;\), (1 < j < il/). 
(4^),- = rf {x? \\M)i){xf >; A,M)P n (xf >; A). 



(3.39) 



They are 'real', S c (x;A) = Bp(x;A) and (x\ A) = Px> jn (a;;A). After some calculation, 
the eigenfunction ( 12.241) is rewritten as 

<j) Vn (x; A) = a I (A' Ml ' MlI ^ , 'a I ^' M '> fn ^ M n K 4«.(^+i)4^n(Mn + i) + lM I M II 

xVto(ar;A)P 2 ,, n (a;;A), (3.40) 

def O (x; A [Ml,Mll] ) [Mj iMn] def 



4>v(x; A) 



Sx,(x - i|; A)Hd(x + i^; A) 



\ + Mj8 l + M U 5 U . 



(3.41) 



The ground state <fix>o is annihilated by Ax>, Ax>(X)4>vo( x ] A) = 0. By construction i/)v(x; A) 
is positive definite in xi < x < x 2 - By using the properties of r)(x), rj(x; A, M), ipm( x ) and 
the determinants, we can show that these Bp (a;) (13.371) and Pr> >n (x) (I3.38P are polynomials 
in the sinusoidal coordinate r)(x): 



E v (x; A) = E v (r](x); X), P Vtn (x;X) = P v , n (r](x); A) 



def 



(3.42) 



and their degrees are generically i and £ + n, respectively (See (IA.6l) - flA.7l) ). Here I is 

Mi Mn 



def 



d f - \ M ^ M i - !) - i M n( M ii - !) + M i M n- 



(3.43) 
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We assume that the parameters are so chosen that Hx>(x; A) 7^ for x\ < x < xi- We have 
also Bp (a;; A + 8) 7^ {x\ < x < x 2 ). By using these and (IA. 16|) f)A. 19[) . we can show that 



3d(x =F %\\ A) 7^ (xi < x < X'l). The lowest degree multi-indexed orthogonal polynomial 
Pr>fi{x] A) is related to Ed(x; A) by the parameter shift A — > A + 5: 

Pv,o{x;X) = AE v {x;X + d), (3.44) 

where the proportionality constant A is given in (lA.ip . This can be shown by using f]A.9[) 
(IA.12I) etc. The potential function Vp (12.28j) after M deletions (s = M) can be expressed 
neatly in terms of the denominator polynomial: 

V v (x; A) = V(x; A^) + - *7i A + 

cx>(a;-i|;A) ^©(x; A + (5) 

The orthogonality relation (I2.23j) is 

/ dxijj v (x; \) 2 P V)n (x] \)Pv t n{x\ A) = /i©,n(A)5 nm (n,m = 0, 1,2, .. .), 

Mi M11 

3=1 i=i 
The shape invariance of the original system is inherited by the deformed systems. The 
operators ^4rf 1 ...d a+1 (A) and ^4d 1 ...d s+1 (A)^ intertwine the two Hamiltonians 7^^...^ (A) and 
"%di...d s+ i(A), 

-4d 1 ...d s+ i(A) t Jl C i 1 ...d s+1 (A) = Hdi...d B (A) - £d s+ i(A), 
A... rfs+1 (A)^ 1 ... ds+1 (A)t = 7^... rfs+1 (A) - 4 +1 (A). (3.47) 

It is important that they have no zero mode, so that the eigenstates of the two Hamiltonians 
are mapped one to one. In other words, the two Hamiltonians Hdi...d 3 (A) and 7^i...d 3+ i(A) 
are exactly isospectral. By the same argument given in §4 of [19], the shape invariance 
of "H(A) is inherited by 7^i(A), 7^^ (A), • • • . Therefore the Hamiltonian 7£d(A) is shape 
invariant: 

A V (\)A V (\) ] = kA v (X + tf)U D (A + 6) + £ X (A). (3.48) 

As a consequence of the shape invariance, the actions of Ap(A) and Aj>(\y on the eigen- 
functions (pv n (x; A) are 

-4d(A)0d„(x; A) = K^/„(A)0xi„_i(x; A + 5), 
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A v (\y(j)vn-i(x; X + S) = k -b n ^(X)(j) Vn (x; A). (3.49) 
The forward and backward shift operators are defined by 

^d(A) = ifo (x; A + 6)- 1 o A V (X) o ^ (x; X) 

= . J , xx + ii; A + *)ei* - S c (a; - A + ty"?*) , (3.50) 

<^(a;)^D(x; A) V / 

Bv(X) = i) V (x; A)" 1 o Av(X) ] o fa (x; X + d) 

V(x; X [Ml ' Mll] )E v (x + q; A)e? p 



E!x)(x; A + <5) 

- K*(:r; X [M ^ M * ] )E v (x - q- X)e-?A V (x), (3.51) 



and their actions on Px> )n (x; A) are 

Fv{X)Pv,n{^ A) = fn(X)Pv,n-i(x; X + 5), iBp(A)Px>,n-i(^; X + 6) = 6 n _i (A)A,, n (x; A). 

(3.52) 

The similarity transformed Hamiltonian is square root free: 

%>(A) = Mx; A)" 1 o Kp(A) o tfj v ( X ; A) = Bp(A) Jb(A) 

= V(x; X [MM ) Sp(x + '^ A ) / e 7 P _ Z v (x-t T ,X + $) \ 
E v (x-q-,X)\ E v (x;X + S) J 

+ V(o;; AW^) ^ (X "^ ;A) (V- - Sp(* + »7;A + *)\ ? (3 . 53) 
^X)(x + z|;A)V ^x)(.x;A + 5) / 

and the multi-indexed orthogonal polynomials Pv,n{x; A) are its eigenpolynomials: 

U v (X)P v ^x- A) = £ n (A)iV(z; A). (3.54) 

Other intertwining relations are (see (IA.16p - flA.19p ) 

^A dl ... d3+1 (X + S)A du .. ds {X) = A dl ... da+1 (X)A dl ... ds+1 (X), (3.55) 

^A dl ... ds+1 (X)A dl ... ds (><y = A dl ... ds+1 (XyA dl ... ds+1 (X + S), (3.56) 



with the potential function V dl ... ds+1 given in ( I2.20P (with s — > s + 1) 
V dl ... ds+1 (x; A) 



^..Asix + q; A) E dl ,„ ds+1 (x - 17; A) 



"d 1 ...d 8 (a; — if ; A) ..d s+1 (a:; A) 

. ^(A^^V 1 ^; A^) : 4 + i is of type I 
X ^ ^(AWjK'^iiAW) : 4+i is of type II ' l ' J 
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where si and sn are the numbers of the type I and II states in {d±, . . . , d s }, respectively. 

Although we have restricted dj > 1, there is no obstruction for deletion of dj = 0. 
Including the level deletion corresponds to M — 1 virtual states deletion: 



Pv, n (x; A) 



= AP v ,, n (x;\ + 5 l ) 

d\, =0 



V = {d\ - 1, ... , - 1, d? + 1, . . . , 4/„ + 1}, (3.58) 
P Vtn (x;\) = BP- Dl n (x;X+ ~5 l ), 

i?' = {4 + i,..., 4. fi + 1, d? - 1, ... , <4 n _i - ( 3 - 59 ) 

where the proportionality constants A and B are given in (1A.2j) - (1A.3j) . These can be shown 
by using ( 13.211) . ( 1A.9I) . ( 1A.11I) etc. The denominator polynomial behaves similarly due 
to ( 13.441) . This is why we have restricted dj > 1. 

The exceptional Xi Wilson and Askey- Wilson orthogonal polynomials presented in [T2~| 
ITU] correspond to the simplest case M — 1, V — {£} of type I, £ > 1: 

&(a?; A) = A + £5 - S l ), 

~i f -(ai + a 2 + nY l : W 

P «- feA » = W*iA + «-*>*{ (^id-^)- :A W ' (3 ' 60) 

As observed in some multi- indexed Laguerre and Jacobi polynomials [lj, it can happen 
that two systems with different sets T> turn out to be equivalent. Namely the denominator 
polynomials with different sets T> may be proportional to each other. For example, the 
denominator polynomial of k deletions of type I virtual states, D\ = {m , (m + l) 1 , . . . , (m + 
k— l) 1 }, and that of m deletions of type II virtual states, T> 2 = {k 11 , (fc+l) n , . . . , (k+m—1) 11 }, 
are related, 

H 2 5 1 (aj; A + mS ) = AE v , 2 (x; A + k8 ) (k,m>l), (3.61) 
where the proportionality constant A is given in (|A.4|) . From this and (I3.45p . we have 

V Vl (x; A + md U ) = V V2 (x; A + kS 1 ). (3.62) 

Therefore these two systems are equivalent under the shift of parameters. Classification of 
the equivalent classes leading to the same polynomials is a challenging future problem. 

For the cases of type I only (Mi = M, Mu = 0, V = {di, . . . , g?m}), the expressions (13.371) 
and (I3.38P are slightly simplified, 

W 7 [g 1 ,...,gj(x;A) = <p M (x)Z v (x;\), (3.63) 
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4(4 A/+1) ;A) 
4(4 M+1) ;A) 



> 6d M > ^-fnlfa; A) = ^ +1 (a:)Pu ifl (x; A) 



^M(M+1) 



g M (4 M+1) ;A) ri(4 M+1) )^n(4 

(M+l) 



M+l) 



s 



(■''•) ' ; A) v\{x 2 ^)P n (x 



;A) 

(M +1 ). A) 



^d M \ x M+l ' 



\\ _I , (M+lUp / (M+l). A N 
^» y ' JVf+1 V Af+1 / J r)l-^)i//'ii i ^* J 



(3.64) 



where r l -{x) = rAx\ A, M + 1). The cases of type II only (Mi = 0, Mq = M) are similar. 
3.4 Analyticity and Hermit icity 

At the end of this section we comment on the hermiticity of the Hamiltonian %t>- The func- 
tions V(x), <f>o(x) and TZ(x) in § 12.21 correspond to Vr>(x] A), 4>t>o{x\ A) oc ipv(x; \)Pr>fi(x; A) 
and Pz>, n (x; X)/Pv,o(x; A), respectively. So the function G(x) ( 12. 7ft becomes (up to an overall 
constant) 



G(x) 



V(x + q- A [Ml ' Mll] )0 o (x + if; A [Ml ' Mn] ) 2 - 



3n(x; A) 2 



•Px>,n(^ + ih A)Pd,, 



z|;A), (3.65) 



and we have 



G(x) - G*(x) = V(x + z|; A [Ml ' Mll] )0 o (x + A [Mi - Mii] ) 2 x 



V(x) 



1 



Hk(x; A 



x — 



V(x) = Pv, n ( x + X)Pv, m (x - i\\ A) - P v ,n(x - z'§; A) A? 



■in \X ~\~ % 2 , A) . 



3d(x; A) 
-• ^ (3.66) 



From ([HUD and (E3ED, this V0 2 , part V(x + if, \ [Ml ' Mll] )(f> ( x + i ^ \ [Ml ' Mu] ) 2 has no poles in 
the rectangular domain D~ r The function V{x) can be divided by 3d(x; A); V(x)/Ex>(x; A) = 
«c^(x) x (polynomial in r)(x)). Thus the potential singularities of G — G* originate from the 
zeros of 3d (a;; A). The left hand side of the condition (12. 8p vanishes because of G(xi + 
ix) = G*(xi - ix), G(x 2 + ix) = = G*(x 2 - ix) for W and G(xi + ix) = G*(xi - ix), 
G(x 2 + ix) = G*(x 2 — ix) for AW, on the assumption that there is no singularity on the 
integration paths. If 3p(x; A) has no zeros in £) 7 , the function G — G* has no poles in 
and the condition (12. 8p is satisfied. If 3p(x; A) has zeros in D 1 , they appear as complex 
conjugate pairs, a ± i(3 (x\ < a < x 2 , < (3 < ||7|), because of 3p = 3d and 3p(x; A) 7^ 
(xi < x < x 2 ). In order to satisfy the condition (12. 8p . the sum of the residues of G — G* 
should vanish. 

The term V v + V% in H v is 



V v (x-X) + V*(x-X) = V(x-\ [MlMl - 



lv(x + q; A) 3p(x - 27; A + 5) 



3d(x-z|;A) 3d(x;A + <5) 
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This does not cause any obstruction for the hermiticity. The potential singularities in the 
interval x x < x < x 2 are (i) V (x; X [Ml ' Mu] ) and V*(x; X [Ml ' Mu] ) at x = x 1 ,x 2 , (ii) zeros of 
^v{% T ^) * n ^ ne denominators, (iii) zeros of E-p(x; A + 8) in the denominators. For case 
(i), the singularities cancel out as in the original case V(x; A) + V*(x; A). For case (ii), we 
can show E v (x =f i\\ A) ^ (xi < x < x 2 ) by using E dl ,„ ds (x; A) ^ 0, E dl „, ds (x; A + 6) ^ 
(xi < x < X2) and (1A.16l) - (lA.19j) . For case (iii), ( 12 .30 j) and ( I3.57P imply that the denominator 
factor Et>(x; A + 6) disappears, namely Su(x; A + <5) does not give singularities. 

Thus the Hamiltonian "H© is well-defined and hermitian, if Ed{x; A) has no zeros in 
or the residues coming from the zeros of Sx>(x; A) cancel. At present we have no general 
proof of the cancellation nor generic procedures to restrict the parameters so that there will 
be no zeros of E-p(x; A) in the rectangular domain D^. Existence of such parameter ranges 
that Ej)(x; A) has no zeros in D 1 can be verified by numerical calculation for small M. 

4 Summary and Comments 

By following the examples of the multi-indexed Laguerre, Jacobi [1] and (g-)Racah [2] poly- 
nomials, the multi-indexed Wilson (W) and Askey- Wilson (AW) polynomials are constructed 
within the framework of discrete quantum mechanics with pure imaginary shifts j3l[26]. The 
method is, as in the previous cases, multiple Darboux-Crum transformations [33l [3H [36] 
by using the virtual state solutions. The virtual state solutions are derived through certain 
discrete symmetries of the original Wilson and Askey- Wilson Hamiltonians and by definition 
they are not eigenfunctions of the discrete Schrodinger equation. The type I and II virtual 
state solutions are introduced f!3.28|) but they are not related with specific boundary condi- 
tions, in contradistinction with the multi-indexed Laguerre, Jacobi or (g-)Racah cases. Main 
emphasis is on the algebraic structure and the difference equations for the multi-indexed W 
and AW polynomials, (I3.52p . (I3.54p etc., which hold for any parameter range. So far we 
do not have a comprehensive method to determine the parameter ranges which ensure the 
hermiticity of the deformed Hamiltonians and thus the orthogonality of the multi-indexed 
W and AW polynomials. The one-indexed, i.e. T> = {£}, £ > 1, of type I are identical with 
the exceptional W or AW polynomials reported earlier [TJl [T9] . 
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Like the other exceptional polynomials, the multi-indexed W and AW polynomials do 
not satisfy the three term recurrence relations. As in the ordinary Sturm-Liouville problems, 
the multi-indexed orthogonal polynomial Px>, n (y; A) has n zeros in the orthogonality range, 
< y < oo (W) or — 1 < y < 1 (AW) (the oscillation theorem). It is well known that various 
hypergeometric orthogonal polynomials in the Askey scheme are obtained from the Wilson 
and Askey- Wilson polynomials in certain limits. Similarly, from the multi-indexed W and 
AW polynomials presented in the previous section, we can obtain the multi-indexed version of 
various orthogonal polynomials, such as the continuous (dual) Hahn, etc. In that sense, the 
multi-indexed Wilson polynomials are also obtained from the multi-indexed Askey- Wilson 
polynomials. Here we briefly discuss the limits to the multi-indexed Jacobi and Laguerre 
cases. In an appropriate limit the discrete quantum mechanics with pure imaginary shifts 
reduces to the ordinary quantum mechanics |42j. Explicitly the W and the AW systems 
reduce to the Laguerre (L) and the Jacobi (J) systems, respectively in the following way 



W : A = — ,— ,g 1 ,g 2 , l = u) 1 + u) 2 , 9 = 9i + 92-\, 

a\a 2 x l 
AW:q X = {-q h \-q h2 ,q 9 \q a2 ), g = g x + g 2 -\, h = h 1 + h 2 -\, q = e^, x = 2x\ 

(a 1 a 2 a 3 a 4 g- 1 )"% AW x c 2 ^ \U\ U' = (p J ) 2 + + - (g + h) 2 . 

v ' 4 sin x J cos^ x J 

(4.2) 

The ground state wavefunction 4>q(x) (13.101) and the eigenpolynomial P n (x) (13. 8p also reduce 
to those of L and J cases after an appropriate overall rescaling. For the deformed systems 
we have the same correspondence under the same limit. For the AW case, the type I and 
II twists (I3.25P reduce to those of the J cases given in PQ, (g, h) — > (g, 1 — h) and (g, h) — > 
(1 — g, h). For the W case, the type II twist (I3.25P reduces to that of the L case given in [I], 
g — > 1 — g. For the type I of W case, there is subtlety because of negative components of 
t T (A) = (1 — 1 — <7i, g 2 ). For example, in order to obtain the limit of the ground state 
wavefunction, we need certain regularization. The limit of type I becomes unchanging g and 
effectively changing x to ix. This corresponds to the type I of L given in pp. Therefore the 
deformed W and AW systems reduce to the deformed L and J systems in pp. The multi- 
indexed Wilson and Askey- Wilson polynomials (13.381) reduce to the multi-indexed Laguerre 
and Jacobi polynomials given in [1] after an appropriate overall rescaling. 
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When all the parameters eij's are real, we have four other twists, 

t ( 13 )(A) = (1 - A 1; A 2 , 1 - A 3 , A 4 ), t (14) (A) = (1 - A 1; A 2 , A 3 , 1 - A 4 ), 
t( 23 )(A) = (A 1; 1 - A 2 , 1 - A 3 , A 4 ), t^(A) = (A 1; 1 - A 2 , A 3 , 1 - A 4 ), 

because of the permutation symmetry of the a^'s. Algebraically, any one of the six twists 
defines a deformed Hamiltonian. According to the parameter configuration, e.g. a\ < a 2 < 
a 3 < a 4 , etc, the compatibility of any two or more twists and the hermiticity of the resulting 
multi-indexed Hamiltonians would be determined. The detailed analysis of these allowed 
parameter ranges is beyond the scope of the present paper. 

With the present paper, the project of generic construction of multi-indexed orthogonal 
polynomials of a single variable is now complete. It is a real challenge to pursue the possibility 
of constructing multi-indexed orthogonal polynomials of several variables. 
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A Several Formulas 



In Appendix we provide several formulas which are not included in the main text for smooth 
presentation. 

First we give various proportionality constants: 
[ A/ i (. „ i ai i i , \ . rr M n, 

,2MiM i: 

xIIh^" 1 )-"!! .-w*/-"-- 1 ) : AW 



Il^iC-oi - «2 + d) + 1) • nJS(-03 - a 4 + df + 1) : W 



A in (jg31D 



? 2MiMii n?i(«i« 2 ?- 



aia 2 g ^ x ) 



A in fl338|) 



_i)M„+i (ai + a2 + n _ 1) n Mi-i d h 



-Oi 



l,aia 2 g 

x(ar 1 a 2 1 a 3 a 4 )^" [IS 9 



2q n - X ) 



(A.l) 



«2 + a 3 + a 4 + 4 + 1) : W 



AW 
(A.2) 
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B in fl339|) 



' - (a 3 + a 4 + n - 1) Y[fX % df(-a 3 ~a i + a 1 + a 2 + df + 1) : W 
(-l) Ml+M »- 1 (a 3 a 4 g ri - 1 )-^(l - a^a^- 1 ) 



x(a 1 a 2 a 3 : 1 a 4 - 1 )5^n-=i/ flI+J ^^ 



A in (jmp 



[3*] 



■1) 



n; = i(-.7r" j n;ii>3 + a 4 - % - a 2 + 2j) 2fc _ 4j -+i 

Ff m J-j)mr-3 n [im], 

iij=iv ^ 1L- =1 (01 + a 2 



AW 



W 



(A.3) 



Lj=i v™i ■ ~z a 3 — a 4 + 2j) 2m _4 :) '+l 

(-a^ 1 a 2 7 1 a 3 a 4 ) fcm g^ (fc - m)(3fcm - (fc - ,n - 1)(fc - ,n+1)) 

n}=i(l nKVr 1 ^ 1 « 3 a 4 g 2j ; g) 2fc - 4j+ i 



TT m (1 _ qj)m-j [fm], _! _! 2 - . 

Li J= i\ y ; 1 1 (aia 2 a 3 a 4 g^;g) 2m _ 4j+ i 



AW 



(A.4) 



Next we give the coefficients of the highest degree term of the polynomials and Px> )n , 
3d (2;; A) = c^{X)r]{xY + (lower order terms), 

Pd(x; A) = n (A)^(x) £+n + (lower order terms). (A. 5) 

For V = {d\, . . . , d l Mv df, . . . , d^}, they are 



Mi, 



c£(A) 



n^w A ))-n^(^( A )) 



n (4-4)- n (4-4 1 ) 

i<j<fc<Mn 



x < 



l<j<k<Mi 

x n*i n£(-«3 - «4 - 4 + ai + « 2 + 4 1 

11.L i /■■ ■ 1 1 — / / 

l<j<fc<Atfn 



n kiK-4)(i_ ? 4-^). n 



i<i<fc<iWi 

lj=l llfc=l ^01020304 



d n -d n \ 



a 3 a 4 g 3 — a x a 2 g 



4,^ 



W 



AW 
(A.6) 



c^,n( A ) = c|(A)c n (A) 



n^i(-«i - «2 - n + 4 + 1) • n*S(-«3 - a 4 - n + 4 1 + 1) : W 
x < g 2M,Af n n j A £ I 1 (aia 2 )"5g^4+ 1 -«)(l - a^^- 1 ) 

x EI •2(«3a4)-^^ ( ^ I+1 ^ ) (l " asa^"- 1 ) : AW 



(A.7) 



The virtual state wavefunction 0x> v (I2.24p for P = {d\, . . . , c?^, ci* 1 , . . . , <i MlI } is given by 

Hx>v(z; A) 



■>Vv 



x; A) 



Sx,(x - A)Hk(x + A) 
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x (a 1 a 2 )^ Ml («3a4)^ /lI K |M(M+1) ^ Ml(Ml+1) ^ A/n(MlI+1) 

^MuiMu-V-M^u^. A W,M„] ) ; v i s f type I 
K kMim-V-MiMnftl( x . A [Mi,M„] ) . v is f type n 

The polynomial £ v with parameters A and A + 8 are related in the following way 
• (vKx; A + ~S 1 )e^ - Vl (x; A + ?) e -^)g(x; A) = # v (A)£(a:; A + 8), 



(A.8) 



ip(x) 



v 2 (x; \)e& - v*(x; A^ A + 6) = %J\)? y (x; A), 



(p(x) 



v 2 [x; 



if(x) 



X + r)ef - ^ 2 (x; A + r) e -**)£(x; A) = /f v (A)^(rr; A + <5) 



«!(x; A)e** - vl(x; A) e -^)ff (rr; A + <5) = ^ v (A)^(x; A), 



(A.9) 
(A.10) 

(A.H) 
(A.12) 



where the functions u i and v 2 are 



: W 



e" ix n j= i(l -%e") : AW' " 



def J nj=s( a i + ia : W 

e-^nLsl 1 -%e") :AW 



(A.13) 



and the constants / SjV and 6 S>V are 
/L(A) = 



a l + 0-2 — v ~ 1 : W 

-q^ v - s \l - a l a 2 q~ w ~ 1 ) : AW 



«3 + 04 — v — 1 : W 

_ g i(v-»)(! _ a 3 a 4 g~ v " 1 ) : AW 



a 3 + a 4 + v : W 

-g-5( v -*)(l-a 3 a 4 g v ) :AW ' 

(A.14) 

ai + a 2 + v : W 

-q-^ Y ~ s \l - a x a 2 q v ) : AW ' 

(A.15) 



These relations are generalized to the denominator polynomials ^di...d 3 v 111 t ne res t °f 
Appendix we consider the set T> = {d\, . . . ,d s }, in which the number of type I virtual states 
is Si and that of type II is sn, s = sj + sn- When v is of type I, the denominator polynomials 
3p v with A and A + 8 are related as 

vl(x; A [SI ' SII] + ~8 l )^ v (x + q-,\ + 8)e^ p 



ip(x)E v (x; A) 



ip(x)E v (x; X + 8) 



- v x (x; X [s ^ ] + S l )E v (x - q- X+ 8) e -'^E Vv (x; A) 
= K-^fl v (X)E Vv (x;X + S), 

v 2 {x; X [s ^ ] )E v {x + q- X)e* p 

- v* 2 (x; X [si ^ ] )E v (x - q; A) e -^)s Cv (x; X + 8) 



(A.16) 
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= K s Hl tV (\)E Vv (x;\). (A.17) 

When v is of type II, they are 

, rr (v* 2 (x; X^ + ~S U )E v (x + if; A + S)e*> 

(p{x)^ v {x; A) V 

- v 2 (x; \ [si ' Sll] + ~S U )E v (x -if,X + d)e-^E Vv (x; A) 

= ^fJ Y (X)E Vv (x;X + 8), (A.18) 

TT^fr , j: (vi(x; X^)E v (x + q- A)e^ 
(p(x)^ v (x; X + 6) V 

- ^(x; A^)H^(x - X)e~^E Vv {x; X + d) 

= kH 1 I v (X)E Vv (x;X). (A.19) 

These relations are used to show fl3.55l) - fl3.56l) . and imply the difference equations of Ex> v ; 
when v is of type I, they are 

(a\X^)V'\x; X^) Ev{X + ^ A) e^ 
V 5j3(z-«2;A) 

+ a\X [s ^)V' u (x; A^) ggfezlil^ e -TP 

^©(x + z|; A) 

E v (x-i^;X) E v (x;X + 6) 
- V*(x; ^u ^-^m^ + irA + Sn 



Er>(x + i%;\) E v (x;X + 8) 
£ l v (X)E Vv {x;X), (A.20) 



and when v is of type II, they are 



e -7P 



a ii (A [si, SII ] ) y> ii (x; X [stM ) Sp(x + ^' A ) e 7P 

- z^; A) 

+ ^(A^^V 11 *^; A^) Sp(X ~^ ;A) 

+ i\\ A) 

- A^" 1 ) ^ X + ^ A) ~ ' 7; A + <5) 
Hx»(a; - £2; A) E v (x;X + S) 

^©(^ + ij, X) ^x)(^;A + o) / 
^(A^faA). (A.21) 
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